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Calculations of two-particle correlations usually assume particles interact only pair- wise after their 
final collisions with third bodies. By considering classical trajectories, we show that interactions 
with the mean field can alter the spatial dimensions of the outgoing phase-space-density profiles by 
tens of percent, consistent with more complicated quantum complications. 
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Two-particle correlations provide insight into the space-time development of relativistic heavy-ion collisions. Speci- 
fially, the six-dimensional correlation function C(P,Q), measured as a function of the total and relative momenta, 
provides information about the the shape of the outgoing phase space density distribution for particles with momen- 
tum P/2. In addition to the size and lifetime of the source, the shape of the phase space packet is also affected by the 
ly-j ■ mean-field traversed by the particles after their last collision. Recently, these distortions were calculated assuming 
T-H [ particles passed through a time-independent optical potential Ji, i^. If the pions left a region of lower mass, it was 
found that the lensing effect of the potential was to distort the extracted source dimensions by a few tens of percent. 
' The effect was analagous to lensing effects from Coulomb mean fields, which were shown to explain different apparent 
^ . source sizes for positive and negative pions at AGS energies 3, 4J. 

The refractive corrections studied in Reference 0, Q allowed for a more physically interpretation of the correlation 
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, data. When ignoring the corrections it appears that the RHIC fireball grows to a transverse (to the beam) radius of 
QQ ■ Ki 13 fm, expanding at a speed « 0.7c, and rapidly disintegrates at a time of « 10 fm/c 0. The puzzling aspect of 
' this picture is that the fireball surface must expand at a speed of 0.7c from the initial collision to grow from its initial 
lO I size of 6 fm to 13 fm in 10 fm/c, not allowing any time for the matter to accelerate. Correlation analyses typically 
provide three dimensions, -Riong) the longitudinal size defined parallel to the beam. Rout, the dimension of the phase 
. space packet perpendicular to the beam and outward along the direction of the pair's momentum, and i?sidc, the 
sideward dimension which is perpendicular to both the pair momentum and the beam axis. The apparent sideward 
dimension in was shown to be increased by the refractive effects of the potential, which would suggest the true 
radius of the fireball was somewhat smaller than the 13 fm previously believed. Reducing the size by one or two fm 



(— I would provide a more physically plausible picture of the reaction's evolution. 
!• • The calculations of Ref. [l|, Ql involved solving for single-particle outgoing wave functions in the presence of 
. 5^ ' complex optical potentials. The wave functions are then symmetrized and the interference term provides the leverage 
. for relating the spatial characteristics of the source to the measured correlation function. Quantum calculations such 
' as these have some drawbacks. First, they are somewhat difficult to interpret due to the inherent complexity of 
solving for wave functions. Secondly, implementation becomes complictated if one were to account for time-varying 
potentials, or potentials without the spherical and boost symmetries assumed in 0,0 ■ Finally, it is difficult to extend 
such calculations beyond the case of identical particles as one then needs to consider the quantum three-body problem. 
For these reasons, it is important to understand the validity which these distortions can be calculated classically. Given 
that the inferred diameters of the RHIC fireballs are often near 25 fm, one might expect that classical considerations 
could be valid except at very low pt . 

In the next section we review the relation between the measured correlation function C(P, Q), the emission proba- 
bility s{p, x) and the asymptotic phase space density profile, /(p, r, i — > oo). By utilizing the connection between the 
phase-space density and the measured correlation, it is shown how classical calculations of the single-particle phase 
space density, which can be readily adapted to include mean-field effects, can be used to generate correlations. To 
provide a physical explanation of how mean-field effects alter the final phase space density and therefore the corre- 
lation function, a simple analytic model of a static cylindrical source is presented in Sec. IIIII Section llVI presents a 
direct comparison of a quantum calculation with one based on classical trajectories. Collective flow plays a critical 
role in correlation phenomenology. A more sophisticated model incorporating longitudinal and radial flow is analyzed 
in Sec. This paper focuses on the distortion of the sideward dimension, which is mainly due to refraction, but the 
apparent outward dimension is also affected. Various issues concerning the outward dimension are briefly discussed 
in Sec. IVII and conclusions are present in Sec. IV 111 
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II. THEORY BACKGROUND 



For non-interacting identical particles, two-particle probabilities can be calculated in terms of one-body source 
functions 6]. Assuming small relative momentum this can be expressed as 
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where the source function s(p, x) represents the probability of emitting a particle of momentum p from space-time 
point X. The correlation due to the interference is the ratio of the first two terms. The source functions can defined 
quantum-mechanically in terms of the quantum T matrices, 



s{p, x) = J2j d*5xT*{x + Sx/2)T,{x - Sx/2)e'P-^'' . 



(2) 



The matrix elements Ti(x) represent the amplitude for a particle to have its last interaction with the source at x while 
the remainder of the source evolves into state i. The last interaction includes interaction through the mean field as 
Eq.s CI and © remain valid in the presence of a mean field, although it is clearly difficult to understand how the 
mean field will alter the T matrices and therefore the source functions. For small relative momentum, one can make 
the smoothness approximation. 



s(P/2, x)S{P/2, y) « s(i?p/2, P/2, a:)s(i?p/2, P/2, y). 



(3) 



which allows one to write the correlation in terms of the outgoing phase space density for particles with momentum 
P/2. 



C(P,Q) = 1 + y"dV5p(r)cos(Q' -r), 



(4) 



5p(r) 



/ d^T'J^r'J{V'/2, , t')f{V'/2, r'„ t')5{v',, - < - r) 
/d3r^d3r^/(P'/2,r;,i')/(P'/2,<,t') 



where the primes denote the positions measured in the pair frame, and Q' = Pa ~ Pj, is the relative momentum in 
that frame. Although S is often called a source function it is more accurately a measure of the outgoing phase space 
distribution for particles of a given momentum. 

An alternative approach is to separate out interaction with the mean field, define the source function so that it 
describes the points at which particles had their last interaction other than through the mean field, then use the 
outgoing wave function, which is a solution to the equations of motion using the mean field, to describe the evolution 
from X to the asymptotic momentum state 0, |3j IM- This involves replacing the phase factor used to describe the 
evolution from x to its asymptotic state p in the derivation of Eq. (Q) , 
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The altered expression is then 
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(6) 



d'^x s{P/2,x)(j)*{pa,x)(j){pb,x) 



where 4> is the outgoing single-particle wave function describing evolution through the mean field. The source functions 
s{p, x) now refer to the points where a particle had its last non-mean-field interaction. 

The motivation for writing these equations is to stress the equivalence of Eq. where the effects of the mean 
field are incorporated into the source function and Eq. ^ where the effects are included by altering the outgoing 
evolution operator. Since both approaches are equivalent, Eq. JB)) can equally motivate the description in terms of 
the outgoing phase space density, Eq. |0J). In this paper we consider applications of Eq. Q where the phase space 
density is calculated by considering classical trajectories through the mean field. 
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ATTRACTIVE REPULSIVE 



FIG. 1: Trajectories that lead to the same asymptotic momentum are altered by an attractive mean field (left). The phase 
space distribution for particles of this given momentum are widened in coordinate space by a factor bmax/R- Repulsive mean 
fields (right) reduce the size of the region of the outgoing particles. 



III. EMISSION FROM A CYLINDER 



As a simple example we consider a pion being emitted from the surface of a cylinder of radius R, where the in- 
medium mass of the pion at the surface is rrimcd and the asymptotic vacuum mass is Wvac- We assume the mass 
returns to its vacuum value exponentially, 

m'ir) = + (mLd " mL)e~^'^-^^/^ (7) 

From time-reversal arguments a thermal source will emit particles of given momentum with the same trajectories 
as those that describe absorption. The asymptotic trajectories of particles with momentum — Pt, Py = 0, that 
intersect with the cylinder have a uniform distribution of impact parameters up to a maximum b as illustrated in 
Fig.Q] The effect of an attractive mean field is to stretch the width of the phase space cloud by a factor 6max/-R- The 
sideward dimension measured in two-particle correlations is stretched by this factor. 

The maximum impact parameter for capture can be found by combining conservation of energy and angular 
momentum, L = bpt, 

pI = pI + K=ff (r) (8) 



2 



Vcnir) = —^+m''{r)-m 

The effective potential will have a maximum outside R unless the momentum is small or if a is large. The condition 
for maximum impact parameter femax is that at the maximum of Vcs defined by 

^ = (9) 
dr 

_ _KnaxPt_ 

the radial momentum is zero so that Eq. (jSJ becomes 



"^vac - m {r)]/a, 



^^*=%^+-'M-™vac- (10) 

Equations Q and can be combined to eliminate and provide a cubic equation for r as a function of b. After 
finding r, one can use either equation to then find pt as a function of b. For higher pt, the location of the minimum 
moves inside R which allows one to ignore Eq. Q and solve Eq. 1)101) with r = R. 

6^^^^i?2!l4ac^!<ed^ P?>(mL-<cd)[W2«)-l]. (11) 
Pt 

The high pt behavior is independent of a, and in the case where a > R/2, the result is independent of a for all pt- 

Figure [5] shows the ratio 6max/ i? as a function of pt for several values of a given the case of a lighter in-medium 
pion mass, rrimcd = 50 MeV/c^. For a = 0, the mean field has no effect as the capture cross section does not extend 
beyond the cylinder. For the saturating value, a = R/2, the cross section is significantly enhanced. For a = R/10, 
the enhancement is identical to the saturating value at large pt and is somewhat reduced at low pt . 
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FIG. 2: Assuming a static cylindrical source, the distortion to the sideward dimension due to the mean field is shown for 
attractive scalar fields. The field lowers the pion mass to 50 MeV/c^ and falls off exponentially outside the emitting radius R 
with a distance scale a. The distortion is stronger for larger ranges a and at lower pt. 



The enhancement of the apparent sideward size is characteristic of an attractive mean field. This is consistent 
with Liouvilles theorem, which states that contraction of the phase space density in momentum space as it leaves the 
attractive mean field should be accompanied by a growth of the phase space cloud in coordinates space. Another 
example of an attractive mean field is the Coulomb field for the residual nucleus which attracts negative pions, but 
repels positive pions. This problem was addressed by Barz, who performed quantum calculations using Coulomb 
wave functions to describe the evolution through the mean field 0, Q . The effect was shown to be significant at AGS 
energies, where it explained the large observed differences between apparent sizes for positive and negative pions. 



IV. DIRECT COMPARISON OF QUANTUM AND CLASSICAL CALCULATIONS 

For potentials where all the characteristic length scales far exceed h/pt, one expects that classical and quantum 
calculations should converge. The radius of the RHIC fireball is close to 10 fm, which readily satisfy the condition 
until Pt is below ^ 50 MeV/c. However, there are other length scales, such as the diffuseness parameter that describes 
how quickly the density falls off a the surface. We consider an optical potential with a Fermi-Dirac form, 

Uir)-{Un + ^UI)^^^::J^J^^, (12) 

where a is the diffuseness parameter. If a = 0, the potential is a step function. Quantum wave packets can reflect 
off step-function potential barriers even when transmission is energetically allowed. This is also true for imaginary 
potential barriers. As an extreme example, a classical particle that heads into a barrier, U{x) = iUjQ{x), with 
Uj oo will be completely absorbed by the barrier, whereas a quantum wave will be completely reflected. Thus, we 
expect the quantum-classical comparison to depend on the sharpness of the potential barrier. 

In order to quantitatively assess the validity of classical trajectory calculations as an alternative to quantum 
calculations, we consider a static two-dimensional complex optical potential as described in Eq. p2|l and solve for 
the distortion to the apparent sideward radii for particles which originate isotropically from points a radius tq from 
the center of the cylinder. For identical particles, radii can be determined from correlations using the expression 



_ l d^C(P,Q) 



(13) 

Q=0 



where the radii represent the dimensions of the asymptotic phase space density, not the dimensions of the emission 
points. Using Eq. (jHl which gives the correlation function in terms of the distributions points for last collisions 
convoluted with outgoing wave functions, one can then apply Eq. (|13|l to write an expression for the variance of the 
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sideward size, 



J Sxs{p, x)<j)* (p, x)<j){p, x) 



where the asyniptotic momentum p moves along the x axis. 

The wave functions (t>{p, x) were calculated numerically by integrating the Klein-Gordon equation for outgoing and 
incoming cylindrical partial waves from r outside the range of the potential to r = 0. Outside the range of the 
potential the partial waves have the form, 

0cii,±(r) -> ^(Mptr/h) ± iYeiptr/h)), (15) 

After integrating the Klein-Gordon equation to find the incoming and outgoing solutions at all r, one makes the 
combination 

Mr) = 'f>iAr)-l3'f>i.-ir), (16) 
where (3 is chosen to satisfy the boundary condition, (f>i{r = 0) = 0. The outgoing plane wave is then 

0(p,r) =0o(r)+ 2i'cosmMr)- (17) 

Depending on the size of ptr/h, the sum over angular momenta might extend to i'max above 50 before results would 
converge. 

Calculations were performed for s being independent of the direction of p with emission points being confined to a 
radius of tq. The potential parameters were = rri^^^^ — m'^^^ with the in- medium mass rrimcd = 50 MeV/c^ and 
Ui = • To with Tq = 100 MeV. This corresponds to a classical decay rate of F = Tom-^/E. The fact that the 
decay rate falls ~ 1/i? is characteristic of a scalar form for the optical potential. The lower panel of Fig. O shows 
the distortion of i?sido, defined as the ratio of i?sidc with the potential to i?sido with [/ = 0, as a function of pt- The 
distance scales for the potential were i? = 10 fm and a — 3 fm. Results are shown for three different values of ro, 5 
fm, 10 fm and 15 fm. The distortions rise at low pt with stronger distortions for small rg. Also displayed in Fig. |2| 
are escape probabilities. Quantum mechanically, the probability that a particle escapes the fireball is 



Jd' ^xsjp, a:)|0(p, r)p 
/ d^xs{p, x) 



-^escape — p ^lA ^~ ( ^ ™\ ■ 



Since low pt particles spend more time in the fireball, and since the decay rate falls as 1/ E, escape rates are small at 
low Pt ■ 

The corresponding classical calculations were performed by solving for the trajectories of particles through the mean 
field. The initial points were chosen on a circle of radius tq with random directions and the energy chosen so that 
their final momentum would have magnitude pt- After the trajectory was calculated, phase space points were rotated 
to make py — Q and {y^) was calculated. Trajectories were calculated consistent with the relativistic equations of 
motion, 

^ = -(l/2i?)V m{r), (19) 

while particles were decayed randomly throughout the trajectory with a rate F = (1/i?) 5t/(r). Classical calculations 
are shown alongside the corresponding quantum calculations in Fig. |3 

The apparent sideward source sizes from the calculations used for Fig. |3| agree well with one another, differing by 
only a percent or two for pt > 100 MeV/c. Even at pt = 40 MeV/c, the calculations agree to better than 10%. As 
explained above, the agreement should be expected to worsen for smaller a. Repeating the calculations for a — 0.5 
fm, discrepancies increased to the level of a few percent for pt > 100 MeV/c, and notably higher for pt ~ 50 MeV/c. 
However, it is hard to physically motivate such a sudden change in densities. 



V. CYLINDRICAL GEOMETRY WITH COLLECTIVE EXPANSION 



The simple cylindrical picture of the previous section ignored collective radial and longitudinal expansion. As 
collective expansion strongly influences both correlation measurements and spectra, it is important to understand 
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FIG. 3: The ratio of the apparent sideward dimension to the dimension without mean field is shown as a function of pt for the 
optical potential described in the text. Distortions are calculated for both quantum calculations (circles) and classical trajectory 
calculations (squares). The source function describing the final collision points was confined to an intial radius ro=5,10 or 15 
fm. The upper panel shows the probability that such particles escape without being absorbed due to the imaginary part of the 
optical potential. The two approaches agree within a few percent. 

how mean-field distortions are affected by both radial and longitudinal expansion. Here, we consider a cyhnder where 
particles are emitted from a surface of radius i? at a time r — \/t'^ — . The matter has a boost-invariant collective 
velocity along the beam axis, — z/t, and a transverse collective rapidity yt at radius R. The break-up surface is 
allowed to move at a different transverse rapidity yf,. This surface should be moving outward, > 0, at early times 
when the system is expanding to reach its maximum radii which are probably near 13 fm. The breakup surface should 
then stop and collapse with becoming negative, and might even become time-like. Thus, the parameters describing 
the expansion, i?, yt, yt and the temperature T, should be chosen based on which phase of the emission one wishes 
to study. 

As in the previous section, we consider a scalar mean field where the pion mass has the form, 

m'{r) = ml, + (m^d " ml,)e-^^-^'+^^*\ (20) 

The initial momenta of the particles are generated randomly to be consistent with a thermal distribution with 
an outward boost of rapidity yt. In addition to the Boltzmann weight, an extra weight proportional to the velocity 
relative to the boundary is included so that the initial generation of phase space points is consistent with the flux of 
a thermal distribution through the boundary. If at any time a trajectory re-enters the interior of the boundary, the 
trajectory is discarded. All trajectories are generated from the same initial point, x — R,y — z — Q,t = t. After the 
trajectory has been calculated, rotational and boost invariance are used to boost and rotate the trajectory so that the 
asymptotic momenta satisfy Py — Pz = 0. The trajectories are then binned as a function of the asymptotic transverse 
momentum, — pt, and the average and variance of the coordinates x,y,z at an asymptotic time t are calculated 
for each pt bin. Since all particles within a given bin have the same velocity, the variances are independent of the 
asymptotic time. 

The variance of the z dimension can be identified as i?ioiig while the variance of y dimension can be identified as i?side- 
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FIG. 4: The distortion of the sideward and longitudinal dimensions due to an attractive scalar mean field are shown as a 
function of transverse momentum. Calculations involved solving classical trajectories through a mean field where the pion mass 
at the breakup surface was 50 MeV/c^ and returned to the vacuum value exponentially relative to the moving surface. Results 
are shown for three values of the exponential scale a. Three sets of parameters are chosen to represent the system while it 
is: expanding (left panels, Vb = 0.4, ii = 9 fm, r = 9 fm/c), has reached its maximum (center panels, Wi, = 0, i? = 12 fm, 
r = 12 fm/c), and is collapsing (right panels, Vh = —0.4, R = Q fm, r = 15 fm/c). For all calculations momenta were initialized 
according to a temperature of 120 MeV and a collective transverse rapidity of 0.8. The distortion is small if the range over 
which the mean field is non-zero is small (a = 0), and becomes very significant when the range approaches 5 fm. Distortions 
are more pronounced when the breakup surface is moving outward as the particles then have more opportunity to interact with 
the mean field. 



Since emission occurs for a variety of radii and times, the experimentally inferred dimensions should be generating by 
averaging over various values for T, i?, t, yt and yi,. For the purposes of gaining insight into the distortions, we study 
the distortions for specific choices of the parameters. The relative weight of emission during the expansion stage, 
where ?/h > and the final burst where j/f, < 0, would depend on details of the equation of state. The variance of the 
X coordinate is also sensitive to the relative times of the emission stages, and we avoid discussing the distortion to 
the outward dimension i?out- 

Trajectories were calculated with and without the mean field. Distortions of the radii are presented in Fig. 01 by 
showing the ratios i?sidc and i?iong to the dimensions calculated without the mean field. For each case the ratios are 
shown for one temperature, T = 120 MeV, and for three values of a: 0, 1 fm and 5 fm. For each choice of parameters, 
one million trajectories were calculated, which was sufficient to suppress statistical fluctuations except at pt below 
100 MeV/c where the pt bins are not well populated. 

The left panels of Fig. 0] show results for an expanding system where both the boundary and the medium have 
positive outward velocities. The surface was assumed to be at a radius of 9 fm at a time r = 9 fm/c while moving 
at a transverse rapidity of j/f, — 0.4. The thermal distribution was characterized by a collective transverse rapidity of 
yt = 0.8. For these parameters, the distortions for i?sidc were even stronger than they were for the static cylinder from 
the previous section. This owes itself to the fact that with yb > 0, the forces from the potential follow the particles 
and push them for a longer time. For this case there is also significant distortion of the longitudinal dimension. The 
center panels illustrate the distortion for particles created when the breakup surface has reached its maximum, yf, = 0, 
with the matter rapidity again chosen to be yt = 0.8. Here the radius and time were both assumed to be 12 fm and 12 
fm/c. The distortion of the longitudinal dimension is small while the distortion of the sideward dimension very much 
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mimics the results for the static cyhnder of Sec. IIIII Finally, the right panels show representative results for when the 
breakup surface is collapsing. Parameters for this case were R — 9 fm, r = 15 fm/c, yt = —0.4 and yt = 0.8. Here, 
the distortions are smaller than for the previous two cases since the particles spend less time under the influence of 
the mean field. 

The results of Fig. 2]can not be directly compared to those of Ref. 0,0 since those calculations used a continuous 
distribution of initial radii modified by an escape probability determined by the imaginary part of an optical potential. 
In their calculation, the in- medium pion mass was also 50 MeV/c^, but since particles tended to be emitted away 
from the center, the mass at the average emission point would be larger. Thus, it is not surprising that the distortions 
of Fig. 2] are somewhat larger, by 50-100%, than those of reference The shape of the pt dependence and the 

relatively larger distortion of i?sidc vs that of i?iong were similar for both calculations. 



VI. THE OUTWARD DIMENSION 



The distortion of the apparent outward size was neglected in this study. The reason for restraining from reporting 
these distortions comes from the nature of the illustrative models applied in the previous sections. In these models 
emission from a specific radius was considered to provide a better understanding of the distortion. Since the outward 
dimension is affected by the duration of the emission, the results for the outward dimensions from these calculations, 
where the radius and emission time are fixed, are unphysical. For instance, if one has a large imaginary potential and 
effectively emits from only the surface, one needs to consistently incorporate the sequential emission from the excited 
regions blocked by the absorptive potential. Thus, any claims about the effect of the real and imaginary potentials on 
i?out can easily be misleading. Nonetheless, we wish to delineate three ways in which the outward size will be affected 
by such potentials. 

(1) A strong imaginary potential confines emission to the surface. Although this might lead to a reduction in the 
-Rout parameter in a toy model, the opposite trend is likely to ensue for consistent dynamical calculations. Pions with 
Pt > 100 MeV/c move faster than the expansion velocity of the fireball's surface. By delaying the emission of pions 
from the interior of the fireball, the pions emitted from the surface get a head start which can lead to the phase space 
cloud being more extended in the outward dimension. 

(2) The particle's velocities are reduced. The group velocity is affected by the medium in both classical and quantum 
calculations. For eikonal approaches such as those discussed in Ref.s 7, 8, 9], the interference between two paths from 
space-time points xi and X2 to asymptotic momentum states pi and p2 is governed by the phase, 

-i{Ei - E2){ti -t2) +i{pi -P2){xi -X2) +i6{pi,xi) ~iS{p2,xi) + S{p2,X2) -S{p2,xi) (21) 
« i{pi -P2) {[{xi - Vpti + (dSi/dp)] - [{x2 - Vpt2 + {dS2/dp)]} , 

where the relative momentum is assumed to be small and the group velocity is dE/dp. Here, S{p,x) is the eikonal 
phase found by integrating over the straight-line trajectory. The derivative of the phase shifts can be equated with 
time and spatial delays, 

/•OO 

d{p,x) = / dx[p{x)-p], (22) 



dS _ 1 
dE ^ 2 



dx 



1 



(23) 



v{x) V 

where p{x) is the momentum of the particle when it passes point x on its way to its asymptotic momentum p. Since 
dx/v is the time spent traversing the interval dx for a classical trajectory, d6/dE is the time delay induced by the 
potential and Vpd6/dE — d6/dp is the spatial offset. Thus, the phase in Eq. (|21|l can be expressed as: 

i{Pi - P2)[xi{t ^ 00) - X2{t ^ 00)], (24) 

where xi — X2 is the asymptotic relative separation of two particles with the same momentum (pi +P2) /2 as calculated 
with classical trajectories. This demonstrates an equivalence (exact in the limit of small relative momentum) agreement 
the eikonal and classical trajectory approaches in a one-dimensional system. The classical perspective provides a 
physically transparent understanding of the distortions to the outward size from a mean field. For instance, the 
apparent outward size can be reduced if those particles emitted early, or further ahead, would be more strongly 
retarded by the mean field than those emitted later. 

(3) Curvature of the trajectories can lead to emission from regions further backward from the direction define by 
the momentum. For the case of an attractive potential illustrated in Fig. some of the trajectories are allowed to 
originate from the edges, and in some cases from the back, of the fireball. For the instantneous-emission pictures 
considered in the previous sections, these distortions tended to extend i?out and i?sido by similar factors. Classical 
Boltzmann equations can consistently incorporate all three of these effects. 
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VII. CONCLUSIONS 

The principal conclusion from this study is that the effects of a particle traversing a mean field between its last 
collision and when it reaches its asymptotic straight-line trajectory can be modeled by calculating classical trajectories, 
then using the final phase space distributions to generate correlation functions. This approach should be valid except 
at very low pt, unless there is a very sharp discontinuity in the mean field. Classical calculations are more tenable 
than quantum distorted- wave approximations. They are simpler to perform and can be applied even when the fields 
are time-dependendent or when there are no symmetries to reduce the dimensionality of the wave-function solutions. 
Another advantage of the classical model is that they can be applied to the case of correlations generated by strong and 
Coulomb interactions by replacing the cosine term in Eq. Qjwhich gives the correlation in terms of the asymptotic 
phase space density, with the squared relative wave function 

The magnitude of the distortions in Fig. 01 are significant, of the order of tens of percent when pt ^ 100 MeV/c. 
However, it should be noted that the conditions chosen for the trajectory calculations in Sec. Ivl were somewhat 
extreme. Not only was the in-medium mass chosen to be have changed by a large amount, the large change was applied 
at the breakup radius from which the particles suffered their last collision. Although this is a strong assumption, it is 
reasonable to assume that the real part of the potential extends beyond the radii at which final collisions occur. This 
is because the real part of the optical potential should be much larger than the imaginary part of the potential at low 
temperature. In RPA calculations, the potentials can be calculated in terms of the phase shifts calculated using the 
relative momenta typical of the medium's temperature. At temperatures below 100 MeV, the invariant-masses of the 
collisions are well below that required for the p and K* resonances, and are beginning to fall below that required for 
the A resonance. The relative strength of the real and imaginary part of the potential is tan (5, and since the phase 
shifts fall well below 90 degrees at low temperature, the real parts dominate. It should be stressed that any attractive 
potential will lead to distortions that broaden the measured sizes. Whereas a falling pion mass is rather controversial, 
in the low density limit one can calculate the mean field with some confidence using measured pion-pion and pion 
nuclcon phase shifts. At very low temperature, only s-wave phase shifts contribute and the interactions are both 
repulsive and attractive depending on the channel. For temperatures in the range of 100 MeV, p-wave interactions 
play an important role. These interactions tend to be strongly attractive due to level repulsion caused by mixing of 
the pion states with A-hole, p-hole or K* states Such mixing was even used to motivate pion condensation at 
high baryon density prj . Though such extreme modifications of the dispersion relation are unlikely given the small 
baryon density at breakup, the overall density of hadrons at breakup is probably near normal nuclear density ~ 0.15 
hadrons/fm'^, and the pion dispersion relation might be modified by several tens of MeV. 

The calculations presented here were only meant to be illustrative. The potential magnitude of mean-field effects and 
the ability of classical pictures to model the effects underscore the importance of incorporating mean field effects into 
hadronic Boltzmann descriptions of the breakup stage. Such calculations have been undertaken in the past, especially 
for lower-energy Fermi- velocity collisions, where the mean field is paramount |l2lll3j| . In this energy range, phase space 
points used to generate correlations are typically taken from when a particle leaves the region of mean field, rather 
than when they have their last collision. Given the likelihood that the dispersion relations for pions are non-trivial, it is 
important that calculations incorporate momentum-dependent interactions. Significant theoretical effort has already 
been applied at developing a theoretical framework for modeling low-density transport that consistently handles both 
collisions and the modification of the particle's group velocities 14, 15, l^[l3- Although this development has largely 
been aimed at intermediate-energy physics, it should be quite applicable for the breakup phase at high-energy. Such 
an effort is crucially important for reconstructing the space-time evolution of the fireball, and for understanding the 
pressure in the hadronic phase. 
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